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The production of a pair of jets with large rapidity separation in hadron-hadron 
collisions, and of forward jets in deep inelastic scattering, can in principle be used 
to test the predictions of the BFKL equation. However in practice kinematic 
constraints lead to a strong suppression of BFKL effects for these processes. This 
is illustrated using a BFKL Monte Carlo approach. 



Perturbative QCD at fixed order in the strong coupling constant a s pro- 
vides sufficient predictive power for a wide variety of high energy phenomena. 
However in some regions of phase space, large logarithms can multiply the cou- 
pling, spoiling the good behaviour of fixed-order expansions. In certain cases 
these large logarithms can be resummed, as in the Balitsky, Fadin, Kuraev and 
Lipatov (BFKL) equation El, and the predictive power of the theory is then in 
principle restored. 

At hadron colliders the BFKL equation applies to dijet production when 
the rapidity separation A of the two jets (with pri ~ PT2 ~ Pt) is large. 
The emission of (real and virtual) gluons in the rapidity interval between 
the two leading jets (see Fig. pi) generates lacge logarithmic contributions 
(a s \n(s/p^)) n ~ (a s A)™ which when resummedel give a dijet subprocess cross 
section which increases with A, ajj ~ exp(AA), with A = 121n2a s /7r rs 0.5, in 
contrast to the &jj — > constant behaviour expected at lowest order. In practice 

"Presented by W.J. Stirling at the 6th International Workshop on Deep Inelastic Scattering 
and QCD (DIS98), Brussels, 4-8 April 1998. 



1 



PT2>+A 



k T ,Xj 



Figure 1: Schematic representation of (a) dijet production with large rapidity separation A 
in hadron— hadron collisions, and (b) forward jet production in deep inelastic scattering. Note 
that the solid lines at each end of the gluon ladder in (a) and at the bottom of the ladder in 
(b) represent either quarks or gluons, according to the effective subprocess approximation. 



such behaviour can be difficult to observe because a gets folded in with parton 
distribution functions (pdfs), which decrease with A much more rapidly than 
the subprocess cross section increases: 

<Tjj ~ q(xi) q(x 2 ) exp(AA) (f) 

with xi } 2 ~ PT/v / sexp(A). The challenge is to find measurable quantities in 
dijet production that are insensitive to the pdfs, but that retain the_distinctive 
behaviour characteristic of BFKL resummation. One possibility uuB is the 
azimuthal decorrelation of the two jets: the multiple emission of soft gluons 
between the leading jets predicted by BFKL leads to a stronger decorrelation 
than does fixed-order QCD, and the prediction is relatively insensitive to the 
pdfs. Another possibility ijLto look for the increase in a with A by considering 
different collider energies □. The idea is to choose A's that correspond to 
the same parton momentum fractions at different energies so that the pdf 
dependence is the same for both, thereby allowing the A dependence of a to 
be extracted. This has been studied recently in Ref. see alsoQ. 

Exactly the same (BFKL) physics applies to the production of 'forward' 
fcy ~ Q 2 jets in deep inelastic scatteringtH, where the large rapidity separation 
is between the current and forward jets, see Fig. [j]b. 

However, a key issue is whether the kinematic regions currently accessible 
at the Tevatron and at HERA allow the asymptotic BFKL conditions to be 
fulfilled. In particular, the asymptotic exp(AA) behaviour is obtained under 
the assumption that there is no kinematic penalty for emitting arbitrary large 
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numbers of gluons with kxi ~ pt- In practice, there is an overall constraint 
from < \fs and even before this is reached there is a strong suppression of 
large subprocess energies from the pdfs. To investigate these effects, we have 
constructed^ a BFKL Monte Carlo programme (BFKL-MC) which reproduces 
the analytic predictions by explicitly taking into account the emission of real 
and virtual gluons to all orders.^] With such an event generator the effects of 
kinematic constraints and experimental cuts can readily be taken into account. 

Fig. H shows the dijet cross section as a function of A for the naive BFKL 
and improved BFKL-MC cases at y/s = 630,1800 GeV, with pti,PT2 > 
20 GeV. Asymptotic QCD LO is also shown for reference. The naive BFKL 
cross section (dashed curve) is always largest, because it includes the analytic 
subprocess cross section &jj ~ exp(AA), which corresponds to the emission 
of any number of gluons with arbitrarily large energies. The prediction falls 
off rather than increases because <x is multiplied by the pdfs, but even those 
incorporate only lowest order kinematics in this case. When exact kinematics 
for entire events are included in both the subprocess cross section and the pdfs, 
as in the BFKL-MC (solid curve), there is a dramatic suppression of the total 
cross section.| In fact the suppression is so strong that it drives the BFKL-MC 
cross section below that for asymptotic QCD LO. The reason is due to simple 
kinematics: the QCD LO cross section contains only two jets, but the BFKL 
MC cross section also includes additional jets, each of which increases the 
subprocess centre-of-mass energy and elicits a corresponding price in parton 
densities. In the naive BFKL calculation, the contribution to the subprocess 
energy from additional jets is ignored and their net effect is to combine with 
the virtual gluons to increase the subprocess cross section. 

Finally we consider forward jet production in deep inelastic scattering at 
HERA, Fig. ^]b. It is relatively straightforward to adapt the dijet formalism 
to calculate the cross section for the production of a forward jet with a given 
Ut and longitudinal momentum fraction Xj 3> a?Bj- In fact there is a direct 
correspondence between the variables: px2 &t and A <-> hi^Xj/xsj)- In the 
DIS case the variable pxi corresponds to the transverse momentum of the qq 
pair in the upper 'quark box' part of the diagram. In practice this variable 
is integrated with the off-shell 7*3* — * qq amplitude such that p\ x ~ Q 2 . As 
a result, it is appropriate to consider values of fc|> of the same order, and to 
consider the (formal) kinematic limit Xj/xBj —> 00, Q 2 fixed. In this limit we 
obtain the 'naive BFKL' prediction oj et ~ (xj/xBj) , the analogue of &jj ~ 
exp(AA). 

6 The technical details can be found in Ref. |. A similar technique has been used in Ref. 0. 
c The running of a s , which we include, also contributes to the suppression, but it has a much 
smaller effect than kinematics. 
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Figure 2: The dependence of the BFKL and asymptotic QCD leading— order dijet cross 
sections on the dijet rapidity separation. The three curves at each collider energy use: (i) 
'improved' BFKL— MC (solid lines), (ii) 'naive' BFKL (dashed lines), and (iii) the asymptotic 
(A 2> 1) form of QCD leading order (dotted lines). 



Fig. H shows the differential structure function d 2 F^j 'dxjdkj, as a function 
of XBj at HERA, with 

Xj=0.1, Q 2 = 50GeV 2 , Q 2 /2 < k 2 T < 4Q 2 . (2) 

The lower dashed curve is the leading-order prediction from the process j*Q — > 
qqQ, with Q — g,q, with no overall energy-momentum constraints. This is the 
analogue of the ajj — > constant prediction for dijet production. Note that here 
the parton distribution function at the lower end of the ladder is evaluated 
at x — Xj, independent of xbj. In practice, when xej is not small we have 
x > Xj and the cross section is suppressed, as indicated by the lower solid 
curve m Fig. 1 The upper dashed curve is the asymptotic BFKL prediction 
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with the characteristic (xj/xBj) x behaviour. Finally the upper solid line is 
the prediction of the full BFKL Monte Carlo, including kinematic constraints 
and pdf dependence. There is evidently a significant suppression of the cross 
section. We emphasise that Fig. || corresponds to 'illustrative' cuts and should 
not be directly compared to the experimental data. Nevertheless, the BFKL- 
MC prediction^ do appear to follow the general trend of the HI and ZEUS 
measurements E9. A more complete study, including realistic experimental cuts 
and an assessment of the uncertainty in the theoretical predictions, is under 
way and will be reported elsewhere 
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Figure 3: Differential structure function for forward jet production in ep collisions at HERA. 
The curves are described in the text. 
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